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The spin polarization versus temperature at or near a fully filled lowest Landau level is explored 
for finite-size systems in a periodic rectangular geometry. Our results at v — 1 which also include 
the finite-thickness correction are in good agreement with the experimental results. We also find 
that the interacting electron system results are in complete agreement with the results of the sigma 
model, i.e., skyrmions on a torus have a topological charge of Q > 2 and the Q — 1 solution is like 
a single spin-flip excitation. Our results therefore provide direct evidence for the skyrmionic nature 
of the excitations at this filling factor. 
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At the Landau level filling factor v = 1 [v = N e /N s 
where N e is the electron number and N s — AeB/hc = 
A/2it£q is the Landau level degeneracy and lo is the mag- 
netic length) the ground state is fully spin polarized with 
total spin S = N e /2 0. In recent Knight-shift spin 
polarization measurements a precipitous fall in the 
spin polarization was observed when either one moves 
slightly away from v = 1 or the temperature is increased 
at v = 1. This effect has also been observed in subse- 
quent experiments with tilted magnetic field as well as 
optical absorption studies S. Theoretically, such a re- 
sult is explained as due to the fact that the low-energy 
charged excitations are spin textures (skyrmions) 
instead of the single spin-flip excitations ffl (the latter 
excitations are, of course, possible only for large values 
of the Zeeman energy, or large values of the ^-factor). 
In fact, finite-size system calculations in a spherical ge- 
ometry indicated that, when one adds or removes a flux 
quantum at v = 1 the total spin changes to S = 

Theoretical studies at v = 1 indicated that Q for large 
values of g the excitations are of single-particle type, i.e., 
they carry charge ±e and spin S z = \ and they have the 
size of magnetic length £ . As g is decreased, the excita- 
tions still carry charge ±e, but they cover an extended 
region and have a nontrivial spin order: at the boundary 
of the system the local spin takes the value of the ground 
state and reversed at the center of the skyrmion. Along 
any radius, the spin gradually twists between these two 
limits. The size of the skyrmion is determined by the 
competition between the interaction energy and the Zee- 
man energy. The former favors a large size in order to 
have uniform charge density, while the latter, with in- 
creasing strength, tends to reduce the size. 



Earlier theoretical studies of spin polarization versus 
temperature at v = 1 involved a continuum quantum 
field theory of a ferromagnet as a model and its prop- 
erties at finite temperatures Q. The other work was 
based on a many-body perturbation theory [Q . A quali- 
tative agreement between the calculated temperature de- 
pendence of the spin polarization from these theories and 
the observed results was achieved. 

We have employed finite-size systems in a periodic rect- 
angular geometry to study the spin polarization (S z ) as 
a function of temperature at and about v = 1. Here 



(S Z (T)) ee i(0|S,|0) 



^ z 



,-e s /kT 
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where |0) is the ground state, Z = J^j e~ l3ej is the parti- 
tion function and the summation is over all excited states 
\j) with energy Ej. The ground state and the excited 
states are calculated from the exact diagonalization of a 
few-electron system Hamiltonian in a periodic rectangu- 
lar geometry ||. 

The results for (S Z (T)) versus T (in units of e 2 /e£o 
where e is the background dielectric constant) in an eight- 
electron system at v = 1 is plotted in Fig. 1(a) where the 
magnetic field is held fixed (B = 10T) but the <?-factor 
is varied (0.1 — 0.5). In recent experiments j9) it was 
shown that applying hydrostatic pressure on the electron 
system, one can vary the g-factor at a given magnetic 
field. The experimental results by Barrett et al. are 
also plotted in Fig. 1 for comparison. The observed data 
show a much sharper drop with increasing temperature 
than the theoretical results obtained here. As mentioned 
above, the size of spin-excitations is dictated by the com- 
petition between the Zeeman energy (which is controlled 
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here by the g-factor) and the interaction energy. Inter- 
estingly, the interaction potential can also be modified by 
including the finite-thickness correction to the electron- 
electron interaction, in the calculation 1 10 , Tlf] ■ This is 
shown in Fig. 1(b) where we present the results for a 
finite-thickness correction parameter (3 = 0.5 as an exam- 
ple (H]|ll|] . The agreement with the experimental results 
now improves noticeably. The changes in our results are 
most pronounced for large g, which is a direct evidence 
of the competition between the Coulomb and Zeeman 
energies mentioned above. Our results are also in good 
agreement with recent magnetoabsorption spectroscopy 
results Jl| (Fig. 1). 

The skyrmion description of quasiparticle (hole) exci- 
tations near v = 1 jl 13 assumes that the low-energy, 
long-wavelength effective Hamiltonian of the system is 
given by 



TL = J dx 7 diii(x) ■ diia(x) 



dx / dy q(x) U(x,y) q(y). 



(1) 



Here n(x) is a unit vector field representing the local 
spin polarization, U (x, y) is the inter-electron interaction 
potential and 7 is a constant which is related to the spin 
stiffness. In our following arguments we shall consider 
the g = case. The deviation of the electron charge 
density from its v = 1 value is given in terms of the spin 
density 



q{x) 



47T 



: n(x) • (diii(x) x djn(x)) 



(2) 



where the right hand side is the topological charge den- 
sity. It is to be noted that Q = J q(x) is always an 
integer which is equal to the number of times n(x) wraps 
around the sphere as x varies over all space. 

In Eq. ([!]), both the terms are small if n(x) is slowly 
varying i.e., if its derivatives are small. We therefore ex- 
pect that smoother the configuration, smaller its energy 
would be. The smoothest spin configuration that one can 
imagine on any closed surface (e.g., a sphere or a torus) 
is the hedgehog configuration. That is the spin configura- 
tion where the spin is always normal to the surface. As 
we go over all space of a sphere the spin density vector n 
covers the sphere exactly once. The charge of this state 
is therefore Q = 1 and the total spin is zero. This effect 
has been observed numerically in the interacting electron 
systems [||. 

Let us consider the electron system on a plane with pe- 
riodic boundary conditions, neglect the interaction term 
initially and consider the energy functional of the sigma 
model [the first term of Eq. ([!])]. All the multisoliton 
solutions of the sigma model are exactly known In 
terms of the variables w — cot(|)e l<6 (which correspond 
to the stereographic projection of the sphere onto the 



complex plane) the solutions are analytic and antiana- 
lytic functions. They saturate the Bogomol'nyi bound 
for the energy which is given by 



E>torf\Q\. 



(3) 



For periodic boundary conditions the solutions are dou- 
bly periodic analytic functions, namely elliptic functions. 
A basic property of elliptic functions is that the sum of 
its residues at its poles in the fundamental domain is 
zero jl5| . The charge Q of the configuration is therefore 
at least two. The winding number is calculated as |l4| ] 
Q = rii, where i runs over all the poles and is the 
order of the i-th pole. Now, if the sum of the residues 
is to vanish, then the elliptic function must have at least 
two simple poles or a second-order pole. Therefore we 
must have Q > 2 for the solutions in this geometry. Pic- 
torially, this would mean that the hedgehog on a torus 
has Q = - the sphere is covered twice as we go over a 
torus, but it is covered in the opposite sense such that the 
total winding number is zero. In other words, the hedge- 
hog on a torus can be viewed as a skyrmion- antiskyrmion 
pair. But if we reflect all the spins in the inner half of 
the torus about the xz (or yz) plane, then the sign of the 
topological charge density flips in the inner half and we 
get Q = 2. The total spin of both the above configura- 
tions is zero by symmetry. 

In what follows, we show that the Q — 1 configu- 
rations that saturate the bound in Eq. (^|) are singu- 
lar and pointlike. Consider the configuration w(z, z) — 
A sn(z) e in ( z ' z \ where sn(z) is one of the Jacobian el- 
liptic functions with periods 4K,2iK' |l^], and A is a 
constant. If we choose 2K = L x and 2K' = L y , where 
L x and L y are the lengths of the sides of the rectangle, 
then sn(z) is periodic in the y direction and antiperiodic 
in the x direction. If £!(z + L x , z + L x ) = Q(z, z) + 7r and 
it is periodic in the y direction then w{z, z) is periodic. It 
is easily verified that if the derivatives of f2 are periodic 
then w(z, z) has Q = 1. The energy of w is 



E[w] = 8irj + 47 / dx 



(A\sn(z)\f 



(l + (A\ S n(z)\y 



In the limit A — > 0, the first factor of the integrand is 
sharply peaked about the pole of sn(z). Therefore, if the 
derivatives of 0, vanish near that point then the bound 
(0) is saturated in the above limit. The Q — 1 config- 
uration that minimizes the energy is therefore singular 
and pointlike. Inclusion of the Coulomb interaction is 
likely to make it regular but it will still be localized in a 
small region near the pole. We therefore conclude that 
the Q — 1 excitation in the rectangular periodic geom- 
etry is more like a single spin-flip excitation and not a 
spin texture. We should point out here that the above 
arguments apply only to the g = case where the solu- 
tions for W are analytic (and antianalytic) and therefore 
sensitive to boundary conditions. 
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We now argue that in the presence of Coulomb repul- 
sion the Q = 2 solutions will have zero spin. We take 
as an ansatz, w(z) = A sn(z), where L x = AK and 
L y = 2K' . The Jacobian elliptic functions have well- 
seperated poles and zeros. Since the poles and zeros 
correspond to north and south poles respectively these 
correspond to slowly varying spin configurations. Fur- 
ther, if we choose A = \fk, where k is the modulus 
p5[ , then using the property sn(z + iK') — [k sn(z)] -1 , 
we have w (z + iL y /2) = [w(z)] , and as a result, 
q(z + iLy/2) = q[z). This choice of A would give the 
lowest value of the coulomb energy since the charge is 
spread out symmetrically. Since cos 9 = {ww — l)/(ww + 
1), the transformation property of w(z) mentioned above 
can be used to show that the z component of the to- 
tal spin vanishes. The x and y components of the total 
spin can be shown to be zero using the property that 
sn(z + 2K) = — sn(z). The periodic rectangular geome- 
try therefore provides a unique test for the geometrical 
and topological aspects of the sigma model scenario. It 
predicts that the spin of the Q = 1 excitations should not 
deviate very much from the ground state spin whereas the 
spin of the Q — 2 excitation should drop to zero. 

Our exact diagonalization studies provide strong sup- 
port for these predictions. In Fig. 2, we present the re- 
sults for (S Z (T)) at v — 1 with one flux quantum added 
(y = 8/9) or removed [y = 8/7). Obviously, the spin po- 
larization does not drop to zero, in contrast to what one 
expects in a spherical geometry, but has the nature of a 
single spin-flip excitation as anticipated above. The rapid 
drop in spin polarization takes place only when we add or 
remove two flux quanta from the system i.e., at v = 8/10 
or at v = 8/6. Indeed, for both fractions (S Z (T - 0)) = 
when g is small (Fig. 3). It should be pointed out how- 
ever, that both v = | and v = | are genuine fractional 
quantum Hall states which were observed experimentally 
jll],[l6| . Additionally, for the latter fraction, the fact that 
it has S = was already established theoretically as well 
as experimentally at T = jy]]. It is quite difficult to 
envision these many-body states as two-skyrmion excita- 
tions of v = 1. Although desirable, it is a formidable task 
to diagonalize even larger systems. Therefore we opted 
for a six-electron system instead. Here, for v = 6/8 or 
v = 6/4 wc do not expect, a priori, that the ground state 
has S = 0. However, as shown in Fig. 4, both these 
states have S = at T = and they also show other 
expected behavior of a spin-singlet state ]s]Ji7p. Obser- 
vation of such a state at these fractions therefore lends 
support to our prediction above that toroidal geometry 
supports at least two skyrmions as lowest-energy spin ex- 
citations at v = 1. 

In conclusion, at v = 1 our finite-size system results are 
in good agreement with the experimental results. Finite- 
thickness correction of the interaction is found to improve 
the agreement, thereby providing strong support for the 



skyrmionic picture of excitations near v = 1. We have 
also demonstrated that the results of the interacting elec- 
tron system are in complete agreement with the sigma 
model predictions. 
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FIG. 1. Electron spin polarization (S Z (T)) as a function 
of temperature T at v — 1 for a 8-electron system without 
(a) and with (b) finite-thickness correction included. Exper- 
imental results of Q (□) and [Q (o, •) are also given for 
comparison. 

FIG. 2. Electron spin polarization versus T at v = 1 for a 
8-electron system with one flux quantum added or removed. 

FIG. 3. Electron spin polarization versus T at v = 1 for a 
8-electron system with two flux quanta added or removed. 



FIG. 4. Same as in Fig. 3 but for a 6-electron system. 
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